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Abstract 

In this paper we study the integrability of the Friedmann equations, when the 
equation of state for the perfect-fluid universe is nonlinear, in the light of the Cheby- 
shev theorem. A series of important, yet not previously touched, problems will be 
worked out which include the generalized Chaplygin gas, two-ternr energy density, 
trinomial Friedmann, Born-Infeld, two-fluid models, and Chern-Simons modified 
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gravity theory models. With the explicit integration, we are able to understand ex¬ 
actly the roles of the physical parameters in various models play in the cosmological 
evolution which may also offer clues to a profound understanding of the problems in 
general settings. For example, in the Chaplygin gas universe, a few integrable cases 
lead us to derive a universal formula for the asymptotic exponential growth rate of 
the scale factor, of an explicit form, whether the Friedmann equation is integrable 
or not, which reveals the coupled roles played by various physical sectors and it is 
seen that, as far as there is a tiny presence of nonlinear matter, conventional linear 
matter makes contribution to the dark matter, which becomes significant near the 
phantom divide line. The Friedmann equations also arise in areas of physics not 
directly related to cosmology. We provide some examples ranging from geometric 
optics and central orbits to soap films and the shape of glaciated valleys to which 
our results may be applied. 

Keywords: Astrophysical fluid dynamics, cosmology with extra dimensions, 
alternatives to inflation, initial conditions and eternal universe, cosmological appli¬ 
cations of theories with extra dimensions, string theory and cosmology. 

PACS numbers: 04.20.Jb, 98.80.Jk 
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1 Introduction 

The Einstein equations of general relativity reduce, when the spacetime metric is of the 
Friedmann-Lemaitre-Robertson-Walker type governing an isotropic and homogeneous 
universe filled with a perfect fluid characterized by its pressure P m and energy density p m , 
to the Friedmann equations: a set of nonlinear ordinary differential equations, determining 
the law of evolution of the spatial scale factor a(t ) , in terms of the Hubble “constant”, H = 
-. It is a challenging task, not always possible, to obtain a complete explicit integration 
in finite terms. In the recent work [T], we demonstrated a link between a theorem of 
Chebyshev and the explicit integration in both cosmological time t and conformal time p 
of the Friedmann equations in all dimensions and with an arbitrary cosmological constant 
A. We were able to establish that for spatially flat universes an explicit integration in t 
may always be carried out, and that, in the non-flat situation and when A is zero and 
the ratio w of P m and p m in the barotropic equation of state of the perfect-fluid universe 
is rational, an explicit integration may be carried out if and only if the dimension n of 
space and w obey some specific relations among an infinite family. Besides, we proved 
that the situation for an explicit integration in p is complementary to that in t so that in 
the flat-universe case with A ^ 0 an explicit integration in p can be carried out if and only 
if w and n obey similar relations among a well-defined family, and that, when A = 0, an 
explicit integration can always be carried out whether the space is flat, closed, or open. 
Furthermore, we illustrated with an example that our method is also useful in the study 
of more realistic cosmological situations when the equation of state is nonlinear. The 
purpose of our current paper to explore the method further to get exact integrations of 


3 



























some other important cases which are not covered in pQ and give a few examples whose 
integrability is beyond the reach of the Chebyshev theorem but can be obtained by other 
means. We shall also see how the integration helps understand some general situations 
when integration is not possible. 

When the equation of state of the perfect-fluid universe is nonlinear, the situations in¬ 
clude the generalized Chaplygin gas, two-term energy density, the trinomial Friedmann, Q 
Born-Infeld, and two-fluid models. Specifically, for the generalized Chaplygin gas model, 
we work on the flat-universe situation with zero cosmological constant and identify all 
integrable cases. For the two-term energy density model, we can do the same. In both 
situations the Chebyshev theorem may be applied directly. In the trinomial Friedmann 
equation model, we show that the Chebyshev theorem is applicable only in the bottom 
three-dimensional space case. For the Born-Infeld type fluid model and assuming a flat 
universe, we show that the Chebyshev theorem works only in the critical coupling situ¬ 
ation when the cosmological constant and Newton constant fulfill a specified condition. 
However, in all non-critical cases, we show that the Friedmann equation allows an inte¬ 
gration for any value of the cosmological constant. Moreover, we study a two-fluid model 
and carry out its integration by the Chebyshev theorem for a closed universe situation. Its 
interest is that in the original setting the integrand is not of a binomial type but it may be 
recast and decomposed into a binomial form so that the Chebyshev theorem is applicable. 
We will also conduct a study of the two-fluid model in terms of the reduced temperature, 
given as the inverse of the scale factor. In this situation the Friedmann equation in con¬ 
formal time does not allow an application of the Chebyshev theorem because it cannot be 
reduced into a binomial form but may still be integrated explicitly when the cosmological 
constant vanishes and universe is either closed or open (the flat case is trivial). The ne¬ 
cessity of introducing several layers of intermediate variables in the process of integration 
often makes it complicated and cumbersome to express the final dependence of the scale 
factor on either cosmic or conformal time transparently, although such a task is always 
manageable. As illustrations, we will work out concrete examples for the models which 
are integrated in terms of cosmic time. These include the Chaplygin gas, Born-Infeld, 
two-term energy, two-fluid, trinomial Friedmann equation, and Chern-Simons modified 
gravity theory models. Although these examples are of different technical features, the 
common lesson gained is that integration allows us to obtain both qualitatively and quan¬ 
titatively accurate knowledge about the solutions, especially regarding the roles played 
by various physical parameters. It should be emphasized that the study on integrable 
cases may often shed profound insight into nonintegrable cases as well. Here we carry out 
such a study for the important Chaplygin gas model. We will see that some integrable 
examples lead us to uncover a universal asymptotic exponential growth rate formula for 
the expansion of the scale factor regardless whether the equation is integrable. The im¬ 
portance of this formula is that it reveals explicitly the roles played by various coupling 
parameters and indicates that, no matter how small the nonlinear Chaplygin component 
is, the conventional linear component always join force to contribute to the exponential 

1 That is, the energy density considered as a function of the scale factor contains three terms, each a 
power of the scale factor. 
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rate, or dark matter, even though the cosmological constant is absent. Moreover, this 
universal formula indicates that, near the so-called phantom divide line, the linear com¬ 
ponent contributes to dark matter significantly. Finally, we will also present a collection 
of examples of the Friedmann type equations which may be integrated in view of the 
Chebyshev theorem but fall out of the subject of relativistic cosmology. 

The plan of the paper is as follows. In Section 2 we introduce the Friedmann equation 
and the relevance of the Chebyshev theorem to its integration. The equation is usually 
expressed in terms of the scale factor a(t), but as, we point out, is sometimes advantageous 
to express the scale factor as a function of cosmic time or to express the so-called reduced 
temperature or inverse of the scale factor in terms of conformal time. In Section 3 we 
review some material about perfect fluids in cosmology. In particular we show how starting 
from an equation of state, we may always express the energy density as a function of the 
scale factor and conversely to obtain the equation of state, given the energy density as a 
function of the scale factor. We also review the Lagrangian approach to perfect fluids and 
its relation to fc-essence and similar models of current interest. In Section 4 we consider in 
detail the Chaplygin gas and the Born-Infeld models. Section 5 is concerned with various 
fluids with multiple components of the type described above. In Section 6 we discuss 
the integration of the reduced temperature form of the Friedmann equation. Finally, in 
Section 7, we point out various occurrences of the Friedmann equation in non-cosmological 
contexts and the relevance of our results to them. In Section 8 we provide a conclusion 
and indicate possible future directions. 

2 Friedmann’s Equation and Chebyshev’s Theorem 

As in pQ, the universe we consider here is an (n + l)-dimensional homogeneous and 
isotropic Lorentzian spacetime with the metric 

ds 2 = g^dx^dx 1 ' = — dt 2 + a 2 {t)gijdx l dx\ i, j = 1,... ,n, (2.1) 

where t is the cosmological (or cosmic) time and gij is the metric of an n- dimensional 
Riemannian manifold M of constant scalar curvature characterized by an indicator, k = 
—1,0,1, so that M is an n-hyperboloid, the flat space M”, or an n-sphere, with the 
respective metric 

g lJ dx\lx J = dr2 + r2 ( 2 - 2 ) 

where r > 0 is the radial variable and dI4 2 _ 1 denotes the canonical metric of the unit 
sphere S 71-1 in M n . Inserting the metric fl2.ljMI2.2p into the Einstein equations 

Gfiu + A g^ v = 8nG n T^, (2-3) 

where G is the Einstein tensor, G n the universal gravitational constant in n spatial 
dimensions, and A the cosmological constant, the speed of light is set to unity, and T fW is 
the energy-momentum tensor of an ideal cosmological fluid given by 

T" = diag {-pm, Pm , ■ ■ ■, P m }, (2.4) 
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with p m and P m the f-dependent matter energy density and pressure, we arrive after some 
standard reduction at the Friedmann equation 


2 \QitG n k 
n{n- 

_ d / 


(2.5) 


where H — ^ is the usual Hubble parameter, / = -gj, and p, P are the effective energy 
density and pressure given by 


p = p H-— P = P _— 

P Pm + 87T G n ’ m 87rG n ’ 

subject to the energy-conservation law 

P + Jl(p + P)R - f>m + + Pm)H - 0. 


( 2 , 6 ) 


(2,7) 


2.1 Reduced Temperature 

It is of interest to consider another variable, called in 
the form 

T=^ 


the “reduced temperature”, of 


( 2 , 8 ) 


where Cq > 0 is a constant, so that at the big bang when a = 0 we have infinite “temper¬ 
ature” and the expansion of the universe leads to a decay of this “temperature”. Strictly 
speaking r is proportional to the redshift 1 + z — and since the wavelength of freely 
propagating particles redshifts by this factor and if those particles are in thermodynamic 
equilibrium then r will be proportional to their temperature T. For the standard case of 
radiation for which the energy density is proportional to a _ ^ r,+1 - ) in n spatial dimensions 
the distinction between the reduced temperature and the temperature of the fluid may 
be ignored. However more generally the two are not the same. For example in three 
spatial dimensions, if P = (7 — 1 )p with 7 a constant we shall see in Subsection 3.1.1 
that p oc Ti 31 while p oc r 37 . Thus only when 7 = |, i.e. for radiation, can the two 
temperatures be identified. 


2.2 Conformal Time 

It is often convenient to introduce conformal time 7 by 

dt = adp. 

We shall adopt the notation f = Then f = ff and (12.51) becomes 

. ,. 2 167rG n 4 . 2 

(a ) = — - -pa — ka . 

n(n — 1 ) 


(2.9) 


( 2 . 10 ) 
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2.3 Chebyshev’s Theorem for Binomial Integrals 

We have seen that integration of (12.51) or (12. 10|) requires knowledge whether the binomial 
integral 

i=S xr[a+m ' (2 - n) 

where p,q,r (r ^ 0) are rational numbers and ct,/3 are nonzero real numbers, may be 
represented in terms of elementary functions. As seen in [I], the Chebyshev theorem 
m si is what is needed to answer such a question which states that the integral (12. lip is 
elementary if and only if at least one of the quantities 


p + 1 
r 


q, 


P+1 

r 


+ q, 


( 2 . 12 ) 


is an integer. 

Our current work, along the line of [T], is to further explore the insight that may 
be achieved by the Chebyshev theorem for the integration of (j2.5[) or (12. 10p and present 
supplementary methods when the problem is out of reach by the theorem. In particular, 
our main focus will be on the situations where the equation of state is nonlinear. 


3 Perfect Fluids 

In this section we recall some standard materials about perfect fluids as used in cosmology 
and their representation, when there exits a velocity potential, in terms of a Lagrangian 
for a scalar field ip with a shift symmetry ip —> ip + constant. 

In n+ 1 spacetime dimensions the energy-momentum tensor is of the form 

= pu^u v + P(g^+ u^u v ) 1 (3.1) 

where p is energy density, P is pressure and tC satisfying g fiu u fl u l/ = — 1 is the (time- 
like) (n + 1) -velocity of the fluid. So as not to clutter up the formula, we include any 
cosmological constant in the total energy density and pressure which, in this section, we 
denote by p and P with the subscript m omitted. For simplicity of exposition, in this 
section, we shall only consider the spacetime dimension n — 3. The extension to all n is 
straightforward. 

3.1 Equations of Motion and Thermodynamics 

Define the projection operator h£ onto the orthogonal complement of the 4-velocity u 11 by 

K = K + => K" = V 1 *" = (3.2) 

Thus for any vector V 

T" = + V £, (3.3) 


7 














with 

vjf* = -u ,t {y v u v ), v» = kv\ 

projecting the equation 

V,,T IW = 0, 

parallel and perpendicular to u M results in 

(p + P) + rPV^p = 0, 


= 


p + P 


hTVvP, 


where 

<P = =>■ = 0, 

is the acceleration of the fluid. The expansion 0 of the fluid is defined as 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


0 = V„u", (3.9) 

and if V is a comoving volume of fluid then 0 = y. 

We now apply the first law of thermodynamics to a volume V of our fluid, assuming 
no conserved charges or chemical potentials. This states that 

d (pV) = Td(sV) - PdV. (3.10) 

Thus 

(dp — Tds)V = —{p + P — Ts)dV. (3.11) 

Assuming that the fluid is homogeneous, i.e., that p and s do not depend upon V, we 
deduce the Gibbs-Duhem relation: 


p + P = Ts. 


(3.12) 


The quantity p + P is the enthalpy density of the fluid and is therefore equal to the 
product of entropy density s and temperature T. Now the first law (13.1111 and Gibbs- 
Duhem relation (13.12p imply 


dP = sdT , <;=> dp = Tds, 


which gives in addition 


dP _ dT 
p + P ~ ~T' 

Thus (13.61) implies the conservation of entropy 

V m (sm m ) = 0, 

and the Relativistic Euler equation (13.81) may now be rewritten as 

= -Grv,;r. 
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(3.13) 

(3.14) 


(3.15) 


(3.16) 











3.1.1 Energy Conditions 

The following three energy conditions are important for cosmology: 

• The Positive Energy Condition: > 0 for every future directed timelike 

vector EC In our context this means that p > 0. Thus for a barotropic equation of 
state such that P = (7 — 1 )p = wp, where 7 is a constant, positive energy requires 
7 > 1 . 

• The Dominant Energy Condition: —TjfV v is a future pointing timelike or null vector. 
In our context this means that p > |P| > 0. For a barotropic equation of state with 
positive pressure, dominant energy requires 7 < 2. The limiting case P = p is called 
stiff matter and this is satisfied by a massless scalar field (see Subsection 13.51) . For 
negative pressure and positive energy, dominant energy requires P > —p or w > —1. 
The limiting case corresponds to a positive cosmological constant 

• The Strong Energy Condition: ( T > 0 for every future directed 
timelike vector CO In our context this means that p + 3 P > 0. For a barotropic 
equation of state, strong energy requires 7 >|orP>—|p. The limiting case is 
attained by an isotropic gas of relativistic strings. Incidentally an isotropic gas of 
relativistic membranes or domain walls has P — —\p and so violates strong energy 
by a considerable margin. 

Note that for all values of k 

- = -^(P + 3P), (3.17) 

a 3 

and so acceleration requires a violation of the strong energy equation. Now 

• 127tG , k 

H = --—(p + P)H—- (3.18) 

O CL 

and so if k< 0 , H > 0 requires p + P < 0. 

3.1.2 Examples 

If P = (7 — l)p with 7 constant then the equation of state is linear and one finds from (13. 141) 

and the Gibbs-Duem relation (13.121) 

s = ir~^T^*. (3.19) 

A simple example of a nonlinear equation of state is provided by the Chaplygin gas 

for which P = — ^ with B constant Mi- One has 

P + P = P--, (3.20) 

P 
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and therefore if B < p 2 , the dominant energy condition, i.e. p > |P|, holds, which is 
consistent with the sound speed v s = ^ being subluminal. The strong energy 

condition, i.e. p + 3 P holds as long as 3 B < p 2 . 

The temperature is given by 

, p 2 >B, (3.21) 

, p 2 < B. (3.22) 

The dependence on the scale factor is 



P 



(3.23) 


where C is a constant. If the scale factor is large, this behaves like dark energy plus (if 
C > 0 ) an admixture of stiff matter. At small scale factors it behaves like dust. 

There exists a generalization [9j [TO] for which 

p = A P ~ (3- 24 ) 

which will be discussed in Subsection 14.11 Note that if a > 0, then at large p the effective 
coefficient 7 = A + 1 and thus A + 1 will determine which if any energy conditions are 
satisfied at large p. 


3.2 Dark Energy 

The simplest interpretation of dark energy is as a cosmological constant term for which 

Tilv = ~8nG 9lJ ' 1 '' (3 ' 25) 

This has perfect fluid form with 


-pDark — PDark — 


A 


8ttG 


= constant. 


Equations f!3.6p and (13.8p are unaffected by the replacements 


(3.26) 


P —■> P + PDark, P —> P — PDark, 9^w 9fiu- (3.27) 
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3.3 Stationary Flows 

For a stationary flow 


y/-g av K"Kv : 

where K IJ ' is a timelike Killing vector field, one has 

u —g av K a K v ), 

whence 

T y/—g av K a K v = constant, 
which is Tolman’s redshifting law. 


(3.28) 


(3.29) 

(3.30) 


3.4 Irrotational Fluids 

The 4-velocity is vorticity free if and only if 


u A du = 0, 

(3.31) 

where u = daT. This implies that 


u = Z~ ?dip 

(3.32) 

for functions Z and ip. Since = 1 we deduce that 


Z = -g^d^ipdvip. 

(3.33) 


3.5 Lagrangian Description of Irrotational Perfect Fluids 

Irrotational perfect fluids admit a Lagrangian description in which the Lagrangian is the 
pressure considered as a function of the velocity potential ip. Such media are often referred 
to as “fc-essence”. The Lagrangian scalar density is given by 

£=y/^L(Z), (3.34) 


whence 


Thus 


2 S/l 

flw = = +2L z d li ipd v ip + g^L. 

p = 2 ZL z — L, P = L. 


(3.35) 

(3.36) 


Because of the shift symmetry ip —> ip + constant, a conserved “enthalpy current” is 


s 


m _ 


1 8C 
\FZ) d tl fl ' 


(3.37) 


Some of examples are: 
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• Stiff matter aka massless scalar L = P = p. 

• Radiation L = Z 2 , P = |p. 

• Polytrope L = Z 2 ^- 1 ') , P — (7 — l)p . 

• Dark Energy L = constant, P = —p. 

• Born-Infeld L = —yT — Z + 1, P = y^. 

• Chaplygin L = —\/l — Z, P = —K 

Note that the last two examples have the same equations of motion but the energy mo¬ 
mentum tensor differs by a dark energy term. That is one gets to the Chaplygin case 
from the Born-Infeld case by the replacements p —> p — 1 and P —>■ P + 1. 

The Polytrope equation of motion is the wave equivalent of the p-Laplacian. That is 

= 0. (3.38) 

If we rewrite (I3.38P as 

d^(^/^{-g ap d a ipdpip) R ^)g ilu d u ip^j =0, (3.39) 

with p = we see that if p — 4, i.e., for radiation. 7 = 1, both the Lagrangian and the 
equation of motion are invariant under conformal rescalings of the metric, i.e., under the 
replacement g )lv —> D 2 g 


3.6 Application to Friedmann—Lemaitre Cosmologies 

In this case d^ip = (ip, 0, 0, 0) and so Z = ip 2 . We regard L(Z) as a function of ip. 
From (I3.36P we have 



L(Z) = P(iP), 

(3.40) 


1 

0, I'i^ 

r 'd 1 

-s>- 

(3.41) 

Now 

^ = (s°, 0,0,0), 

(3.42) 

where 

dip 

(3.43) 

and conservation of entroDV (13.15P 

gives 



,dP 

a —- = constant. 

(3.44) 


dip 


Probably it is best to define say x — ip. 
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One has 


P = P(x), 

ip p 
P = x ^- p 

,d P 

a — = c, 
ax 


(3.45) 

(3.46) 

(3.47) 


where c is a constant. For convenience, we denote g(x) = ^(x) and assume g is invertible. 
Hence we have x = g~ l (ccT z ) and thus 


p = ca 3 g 1 (ca 3 ) 


S(0d£-P(0). 


(3.48) 


This equation relates the energy density p to the scale factor a in a complicated way. 
Nevertheless, many concrete examples are manageable, with or without the Chebyshev 
theorem. Below are a few of these. 


Set P = One finds 

r +1 


p = 


r + 1 


x r+1 = rP. 


(3.49) 


From [I] we know that the Friedmann equation is integrable if and only if - is among 
the sequence of numbers 


- 1 , 


2 5 1 7 4 3 

,_ 3’ _ 9 ,_ 2’ _ 15 ,_ 9 ,_ 7’ 
In particular, the case r = 3, that is 


2 1111 
9’ _ 5’ _ 6’ _ 9’ 3 


ax = constant, 

is well known and corresponds to radiation. 
Chaplygin is 


(3.50) 


(3.51) 


L = -VT 


1 


x x 


P = —, 

P 


and 


x 


Vp 1 


= constant. 


x 


(3.52) 

(3.53) 


A general form of this model will be studied in Section 14.11 
Born-Infeld is 


L = —y/l — x 2 + 1, P — 


P 


P+ 1’ 


and 


x 


vT^ 


= constant. 


x- 


(3.54) 

(3.55) 


The general situation for this model is thoroughly considered in Section 14.31 
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• A model of current interest is the Bose-Einstein condensate considered by Chavanis 

HUS]. 


One has 


P — \/ PoP + P, 

(3.56) 

where Pq is a constant. We have 


dL T 

P = x- -L, 

dx 

(3.57) 

and hence 

dL dx 

2 L + \/PqL x 

(3.58) 

Thus 


L = - i (Ax-P 0 f 

(3.59) 

We choose the constant of integration to agree with the Lagrangian for stiff matter 

so that 


P= l {V2x-P,)\ 

(3.60) 

It follows that 


a?\fP = constant, 

(3.61) 

and hence 

Pdust 0 . Pstiffmatter 0 

P= 3 + 6 

gC a b 

(3.62) 

which is a mixture of pressure matter and stiff matter as stated in 
case of the model investigated in Section 6. 

[4] and a special 


4 The Chaplygin Gas and the Born—Infeld Models 

In this section we study the generalized Chaplygin gas and Born-Infeld models. In the 
former situation, we derive exponentially expanding solutions and obtain their exponen¬ 
tial growth rates explicitly in terms of various physical parameters. In particular, these 
expressions reveal the specific roles played by the linear terms, when nonlinear terms are 
present, in the equations of state. In the latter situation, we show that, although the 
Chebyshev theorem is not applicable in general, an integration can still be performed to 
render explicit solutions. 

4.1 The Chaplygin Gas Model 

The equation of state of a generalized Chaplygin gas model reads 

B , 

Pm = Ap m —, ( 1 - 1 ) 

Pm 
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where A, B are positive parameters and the exponent a stays in the interval 0 < a < 1. 
Inserting (14. lj) into (12.7p and integrating, we obtain 

(1 + A)p" +1 = Ca~ n{ - l+A){ - a+1) + B , (4.2) 

where C > 0 is an integration constant. In fact, in our study below, the positivity 
condition on A can be relaxed to 

A > -1. (4.3) 

So (14.31) will be the condition on the constant A we observe in the subsequent discussion. 
Note that the borderline A = — 1 is known as the Phantom Divide Line (PDL) which is 
an actively pursued topic HD E21 m US □! in the research on dark matter. 

In view of (12. 6 p and (14.20 . we see that (12.5p takes the form 

a 2 = 167rGn (i + yl)-sTTa 2 (C'a- n(1+A)(Q+1) + B)^ + — a 2 - k. (4.4) 
n(n — 1) n[n — 1) 

In the binomial case, A = 0, k = 0, integrability by virtue of Chebyshev’s theorem is 
always ensured. As an illustration, we proceed to carry out the integration. 

For convenience, we rewrite the Friedmann equation as 

a 2 =/3 2 a 2 (Ca~ n{1+A){a+1) + B)^Ti, (3 2 = ^%- (l + A)~&, (4.5) 

n[n — 1) 

whose integration amounts to computing 

/ = J a" 1 (Ca- n(1+A)(Q+1) + B)~W) da. (4.6) 

Since the constants B and C are positive, we may introduce a new variable u such that 


Ca - n (l+A)(a+ 1 ) +B = m 2( Q +1)_ 

(4.7) 

Hence we obtain 

2 f u 2a 

T n( 1 + A)/3 J n 2 ( Q+1 ) - B ’ 

(4.8) 

which may be integrated to yield an elementary function if and only if a 

is rational, say 

mi 1 0 O 

a = —, m-i, m 2 = 1, 2, 3,.... 
rn 2 

(4.9) 

Assuming (jd.Op, we may use 

1 1 

(4.10) 

V — B 2 ( m l+ m 2)'U m 2 


to recast (14.81) into the form 


I 


2rri2 

T n(l + A)(3 B 


m 2 

2 (m^ +7712) 


y2mi+m2 — l 
v 2(mi+m.2) _ \ 


(4.11) 
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and then conduct its integration. 

If a is irrational, we may approach a by a sequence of rational numbers and use the 
above integrals to approximate the desired solution. 

As an illustration, we consider the classical Chaplygin gas model pun] given simply 
by the equation of state 

P m = Ap m - —, (4.12) 

Pm 

which corresponds to m\ = m 2 = 1. Hence (14 . 11 j) gives us the result 


/ 


=F- B 4 

n(l + A)/3 

ip- - - B~* 

n{ 1 + A)/3 



+ arctari v ) + C\, 


(4.13) 


where C\ is an integration constant. The big bang universe a(0) =0 leads to the solution 

2n(l + A)pBk = l J ( B + Ca - M1+A}) ] +B ] 

\(B + Ca~ 2 n ^+A))j - B i 
which may further be resolved to yield 

a 2n{1+A \t) = - 


larctan ^1 + —a 2n ( 1+A )^ + q r, 

(4.14) 


C 


B (cotlr 


where 


n( 1 + A)/3 Bit + b(t) 


b(t ) = arctan ( 1 + ^ a - 2?l ( 1+ ^(t)^ ' - 
B / 2 


(4.15) 


(4.16) 


Since b(t) given in (14.16P is a bounded function, we see that the large-time behavior of 
the scale factor obeys the asymptotic exponential growth law: 


“ W ~ * SB 


(J \ 2n(l + A) 




, t —> OO, 


(4.17) 


which clearly spells out the roles played by various parameters in the model since from 
the above the asymptotic exponential growth rate is seen to be 



nG rj 


n(n — 1) 


1 

2 


B 


1 + A 


1 

4 


(4.18) 


As another illustration, we consider the case when a = \. Thus, inserting rri\ = 
1 , 777.2 — 2 into (14.lip , we have 


4 

/ = T- B 3 

■n(l + A)j3 



dv. 


(4.19) 
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To proceed further, we set w = v 2 * in (14.19ft and use the partial fraction decomposition 


w 


w — 1 


w 3 — 1 3 \ w — 1 w 2 + w + 1J ’ 


to arrive at 

/ = + 


iH fin <” 2 T )2 . + 2V3arctanjd±l'| +Cl , 


3n(l + A)/3 V n 4 + v 2 + 1 1/3 

where C\ is an integration constant. As before, the big bang solution reads 


(4.20) 


(4.21) 


/ r 


3n(l+A)/3B3t = In 


1 + §a~ 2 "( 1+a ) + 1 + %a~ 2 n ( 1+A ) 


\ 


+ 1 


V 


1 + §a"2 n(1+j4) 




1 


—273 arctan —= I 2 

73 


1 + —a~ i n(1+j4) 

B 


+ 11+ T+r. (4.22) 


Resolving this relation, we obtain 


1 + § a ~^ (1+A) w) 3 = 


2 -f e -3n(l+A)l3Bh+2b(t) _|_ ^n{l+A)0B^t+b(t) ^ — 3 e -3n(l+A) / 3BlJt+2b(t) j 2 

2 — e -3n(l+A)/3B3f+26(p 


where 


bit ) = 73 arctan -j= I 2 


C -§n(l+A) 


f + -a-. 




3 \ 73 

+ 1- 71 

' 2 


(4.23) 


(4.24) 


is a bounded function. Therefore we obtain the asymptotic exponential growth law as 
before: 


3n(l + A) 1 

Moo. 


a(t) ~ ( 77 ) 

Thus we also get the asymptotic exponential growth rate as follows: 

7 tG„ V / R \ ^ 

/5Ra =4 


(4.25) 


n(n — 1) J \ 1 + A 


(4.26) 


For later convenience, we use (14.7p and (14.101) to express the dependence of a and v 
directly as 


v = B 


V s ( n -n(l+A) (mi+m2) . 

2(m 1 +m 2 ) I (Jq ' ’ m 2 + iJ 


2(m 1 +m 2 ) 


(4.27) 
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Thus the big bang initial condition a(0) = 0 with A > — 1 leads to 


v(0) = oo. 

Hence, from (14.lift and (14.28ft . we obtain the expanding solution 


n( 1 + A)j3 _ 122 _ 

-1--- ]3 2(m 1 +m 2 ) t — 

2 mo 


fOO v 2mi+m2-l 


’v(t) 


; 2(mi+m 2 ) _ 


dv, t > 0, 


(4.28) 


(4.29) 


which monotonically climbs to the limiting level v = 1 as t —> oo. 

As yet another example, we take mi = 1 and m 2 = 3 (i.e., a = -|). Thus (I4.29P yields 


8 


n(l + A)/3B s t = In 


v + 1 
v — 1 


+ 4 arctan v + \/2 In 


; 2 - + l' 


v 2 + \[2v + 1 

—2v^2 ^arctan + 1 + arctan \/2u — 1 ) 

— 1)7r, n = u(t) >1, t > 0. 


(4.30) 


Since all the terms except the first one on the right-hand side of (14.30(1 stay bounded, we 
have 

v(t) ~ coth ( -n( 1 + A)j3B ® t ] , t —> cxd. (4.31) 


In view of (14.27ft and (14.31ft . we arrive at 


a(t) 


C \ 4 ™(!+ A ) 




t oo. 


16B / 

Consequently, we get the asymptotic exponential growth rate as before: 

3 


(4.32) 


/TBs = 4 


7rG r 


n(n — 1) / \ 1 + A 


B 


(4.33) 


4.2 A Universal Formula for Asymptotic Exponential Growth 
Rate and Cosmological Consequences 

With the explicit knowledge on exact solutions described earlier as stated in (14.1711 . (14.25ft . 
and (14.3011 . we are able to derive such precise formulas as (14.18ft . (14.26ft . and (14.33(1 for 
the exponential growth rate of the scale factor. In these examples, it is interesting, and 
perhaps surprising, to see that, in the context of the Chaplygin gas model, the linear part 
of the equation of state switched on by the coefficient A in (14.ip also contributes to the 
asymptotic exponential growth rate of the scale factor, as seen in the examples (14.18ft . 
(14.26ft . and (I4.33P in the previous subsection, which follow a clear pattern of the following 
characteristics: 

• The growth rate contains a common factor which depends on the Newton constant 
G n and spatial dimensionality n only. 
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• The growth rate is proportional to some power 7 (say) of the common ratio 
composed from the linear and nonlinear coupling constants A and B in (14. ip . In 
(I4.18p . (I4.26h . and (14.33ft . we have 7 = and |, corresponding to a = 1, and 
in (14.ip . respectively. 

Inspired by the above examples, in this subsection, we will derive a universal formula 
for the asymptotic exponential growth rate of the scale factor for the generalized Chaplygin 
model (14.ip . regardless whether a is rational. We then elaborate on some applications of 
this general formula to the study of dark matter. 

For convenience, we use the new variable 


1 

w = B 2 ( q + 1 ) u 


(4.34) 


in (14.8p . Thus (14.7p gives us 

-n(l+A)(a+l) + ^ _ ( 4 . 35 ) 

Therefore the boundary condition consisting of a(0) = 0 and a(oo) = 00 becomes 



tc( 0 ) = 00 , u>(oo) = 1 . 
From (14.81) and (14.361) . we arrive at 


1 1 

-n(l + A)/3B 2 ( a+1 '> t = 


w 2a 


’w(t) 


yj2(a+ 1) _ ^ 


die, w(t) > 1 . 


(4.36) 


(4.37) 


Since u> 2 ( a+1 ) — 1 is approximated by 2 (a + l)(tc — 1 ) for w near 1 , we have 


po ° w 2a 


dw = 


rw(t 0 ) 1 


lw(t) w 2 ( Q+ b - 1 " 2(a + 1) J w(t) w-l 


dw + 0(1), 


(4.38) 


where t > t 0 > 0 and 0(1) denotes a bounded quantity. Combining (14.371) and (I4.38p . we 
obtain 

W(t) = 1 + (W(t 0 ) - l) e ^+A)(a+l)PBm t+ 0(1) ? t > ^ (4>39) 

Using (I4.35P and (I4.39p . we have the asymptotic growth pattern 


i(t) ~ C , 0 e /3s2( “ +1) \ t> t 0l 


(4.40) 


where Co > 0 is some constant. Consequently we see that the asymptotic exponential 
growth rate of the scale factor a is given by the universal formula 


13 B 2 («+i) 


/ 7T G n \ 2 / B \ 
\n(n — 1)/ \1 + Ay 


(4.41) 
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which recovers the formulas (14.18I) . (14.261) . and (j4.33[) when a = 1, and respectively, 
and exactly confirms the characteristics described at the beginning of this subsection. 

ft will be interesting to compare the above result with the situation when the nonlinear 
term in the Chaplygin law is turned off, by setting B — 0 in (14. lj) . Now the conservation 
law (12.71) leads to 

p = Poa - n{1+A \ A e R. (4.42) 

Hence, integrating the Friedmann equation (for A = 0, k = 0), we have the (well-known) 
expanding solution: 


/ 

a(t) = < 




(47rG n p 0 fe] (1 + A) 2 ) HiW i„ ( iU , 

a(0)e 4 V / ^G„po/n(n-l) t, 


( n(l+A) 

a 2 


(0) — 4 


irGnPO 

n{n— 1) 


2 

n(l+A) 


1 


A > —1, t > 0, a(0) = 0; 

A = -1, t > 0, a(0) > 0; ( 443 ) 

A < -1, £ > 0, a(0) > 0. 


So, when the nonlinear matter is absent and linear matter is that of a conventional type, 
A > — 1, the scale factor of a big-bang universe fails to grow exponentially but instead 
follows a power-function type growth law, and, when A < — 1, the big-bang solution with 
a(0) = 0 is not allowed and the universe may only evolve from a positive value of the 
scale factor, a(0) > 0. In this latter situation, the condition A < — 1 corresponds to a 
hypothetical model called the phantom-energy universe PUE] where the solution blows 
up at finite time 


to — 


n(l+A) 

a 2 


( 0 ) 


4 


1 1 

f TTGnpO_\ 2 

y n(n—1) J 


(4.44) 


known as the Big Rip [T8], 19] epoch. The limiting situation, A = —1, is the PDL referred 
to earlier, where the exponential expansion rate of the scale factor may assume arbitrarily 
large value, depending on the value of po- 

On the other hand, the universal formula (14.411) indicates that, within the context of 
the generalized Chaplygin fluid (14.11) . a conventional (A > —1) linear term alone does 
not give rise to exponential growth but it plays a role in exponential growth when a 
nonlinear term is present, no matter how weak this nonlinear matter is, and that, this 
linear contribution becomes significant near the PDL limit. That is, the asymptotic 
exponential expansion rate of the scale factor may assume arbitrarily large values for 
A > — 1 and A stays near —1. In other words, a tiny presence of the Chaplygin matter 
would switch on the presence of dark matter which becomes significant near the PDL 
limit. 


4.3 The Born—Infeld Type Fluid Model 


In this subsection, we shall consider a Born-Infeld type fluid whose equation of state is 
of the form 


P 

1 m 


Pm 

Pm + 1 


(4.45) 
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As we saw in Section [3j there is a close relation between the Born-Infeld and the clas¬ 
sical Chaplygin gas models. The two Lagrangians differ by a constant which adds a 
cosmological or dark energy term. Thus if p m = Pm + 1 and P m = P m — 1 we find 

Pm = ~ (4.46) 

Pm 

However, the Born-Infeld equation of state, unlike the Chaplygin gas, satisfies both the 
dominant and strong energy conditions for all values of the energy density. 

Inserting (j4.45j) into (12.711 gives us the following dependence relation between the mass 
density and the scale factor: 

Pm — — 1 + \/l + Coa ~ 2n , (4.47) 


where Co > 0 is a constant. Hence, in view of (12. 6p . the Friedmann equation (12. 5 p becomes 


H 2 


167rG n 
n(n — 1) 


{—1 + \J 1 + CqCL + 




(4.48) 


First, we see that, in order to be able to utilize the Chebyshev theorem for integration 
of (I4.48p . we have to assume k = 0 and 


A = 8 t rG n . (4.49) 

In such a situation, the Friedmann equation (I4.48P is reduced into (14.51) which corresponds 
to the Chaplygin gas model (14.11) with a — 1, A — 0, B — 1 and whose solution is given in 
(14.13p with C = Co in (14.71) . Thus various solutions of specific interests may be written 
down immediately. For example, the solution satisfying a(0) = 0 is given by 


In 


u — 1 
u + 1 


2 arctan u 





u = (1 + C Q a~ 2n ) * . 


(4.50) 


Next, we note that, when (I4.49P is not fulfilled and, so, the Chebyshev theorem fails 
to be directly applicable, we can still integrate (I4.48P at k — 0. For this purpose, we set 



_ A 2 _ 167rG n 

87 t G n ’ n(n — 1)’ 

(4.51) 

and rewrite (|4.48P as 


(a) 2 = (3 2 a 2 ( a + a/1 + C 0 a 2n ) , 

(4.52) 

whose integration reads 

/ = 

C ' \ i 

/ a -1 [a -\- \/l + C 0 a _2n j da = ±f3t. 

(4.53) 


Obviously we cannot apply the Chebyshev theorem as before. Nevertheless, we may set 
1 + C 0 a~' 2n = u 2 so that a = C 0 2 ™ (u 2 — 1)~^ to convert the left-hand side of (14.531) into 


I 


1 

n 


u(u z — 1) L (a + u) 


'2 d u. 


(4.54) 
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Now set a + u = w 2 or u = w 2 — a. Then (I4.54[) becomes 


/ = 


2 

n 

1 

n 


w 


a 


w 2 — a) 2 — 1 
1 


du> 


+ 


w 2 — (a + 1) w 2 — (a — 1) 


dw, 


(4.55) 


whose integration up to an additive integration constant is readily computed to give us 
the results: 


|o;| < 1, 
a — 1 , 
a > 1 , 
a = — 1 , 

a < — 1 . 


I = 



(4.56) 


Thus, in short, we conclude that the Born-Infeld fluid model (I4.45ji allows an exact 
integration when k = 0 for any value of the cosmological constant, although it lies beyond 
the scope of the Chebyshev theorem. 

To end this section, we examine the case a > 1 closely as an illustrative example. It 
is readily seen that the big bang solution, with a(0) = 0 and a(t) > 0 for t > 0, is allowed 
which is given by 


2 n(3t 


\/l T OL 
1 

y/a — 1 


\J OL yj 1 + Cq(X 2n ' + \/l + OL 

\Jol T \/l T - Cod — y/lTa 

\J ol -j- \/l T CqCl = T \j ol — 1 

L = - 

ol T \/l T Cod ~ !l — v ol —T 


(4.57) 


Thus, in the limit a(t) —> cxd as t —> oo, since the second term on the right-hand side of 
the above is bounded, we have the asymptotic formula: 


a(t) 


G) 


1 

2 n 


8(1 + a) 


7rG n Co \ 2 " 


f)y/T+Ctt 


A 


/ 2A 

n ( n — i) 


t —> oo. 


(4.58) 


It is interesting to note that the Newton constant is absent from the exponent of the 
asymptotic growth law. 
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5 Multiple Fluids 

In this section we study the integrability problem for several multiple fluid models in view 
of the Chebyshev theorem. In a two-term energy density model, we give an example of 
explicit solution that obeys a linear asymptotic growth law. In a two-fluid model, we show 
how to use the information from integration to obtain periodic solutions linking big-bang 
to big crunch epochs. In a trinomial Friedmann equation model, we find exponentially 
expanding solutions. 


5.1 A Two-Term Energy Density Model 

Consider now the classical Friedmann equation in three-spatial dimensions so that (j2.5(1 
assumes the form 


H 2 = 


87 tG 


3 P a 2 ' 


(5.1) 


In 


the energy density is taken to enjoy the general decomposition 


OO OO 

P = £ P>~"+ y p-a-™, 

n=—o o m=—o o 


(5.2) 


where the quantities p+’s are positive constants and, p“’s, negative. Following [20], 2Tj . 
we study the simple situation that only two terms in (15. 2 j) are present and k is set to zero 
so that the Friedmann equation becomes 


H o 


= kl + n a 2 ~ n + fi- a 2 ~ m , 


(5.3) 


where 12+ > 0 and 12“ < 0 are the normalized positive and negative energy density 
coefficients and H 0 > 0 is the Hubble parameter at a time when a = 1, although our study 
here does not need to be confined within such a sign condition regarding the quantities 
12+ and 12“. Hence, from (15.3p . we arrive at the integral 


i = jtJ +S 2 „a”- m r 4 do = ±t. (5.4) 

When m = n, (15.41) is trivially computed. Thus, in the sequel we assume m 7 ^ n. In view 
of the Chebyshev theorem we see that (15.41) may be integrated if and only if m, n (m 7 ^ n ) 
satisfy one of the following conditions: 

2 N 

n = - m, N = 0, ±1, ±2,..., (5.5) 

2N — 1 v 7 

2N 

m = - n, N — 0,±1,±2,.... (5.6) 

2N +1 >11 v ; 


In particular, if m = 0 or n = 0, integrability is ensured for any integer n or m. 
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As an illustration, we show how to integrate (j5.4[) under the condition (15.5I1 . In fact, 
in this case, H5.5[) may be rewritten as 


1 f mN i / , m \ 9 

1 = W 0 J a2N ~ X +^ m a 2JV - 1 ) da. 


Now set 


_ . m o 

+ ^a^-i = u 2 , 


which leads to the inverse relation 


a = 


^ m 


so that (15.71) reduces into the integral 

2(2iV — 1) 


I = 


mHo(n, 


-)N 

m) 


ld “. 


(5.7) 

(5.8) 


(5.9) 


(5.10) 


that allows a recursive computation with the substitution u — . [ sec 0, say. Note that, 
from (15.8p . we have the restriction 

n+-u 2 = \n~\a^ >0 (5.11) 


for the range of the variable u. 

After the above general discussion, we consider a concrete situation with m = 3, n = 
2,N = —1. Thus, from integrating (j5. lQj) . we get 


1^3 I 




where C is an integration constant. On the other hand, since N = — 1, 
strengthen (15.111) into 

— u 2 = |f2j|a -1 > 0. 

In particular, we have the lower bound 


(5.12) 

we need to 

(5.13) 


n(0) > 


1^3 

0+ 


(5.14) 


for the initial value of the scale factor so that a big-bang universe is not allowed. Fix the 
constant C in (15. 12[) with the initial value of the scale factor at t — 0 satisfying (15.141) . 
Choose the plus sign in (j5. 12j) representing an expanding universe. Then, from (j5.12p . we 
see that 


lim u(t) 

t—> OO 



(5.15) 
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or a(t) — * oo as t —* oo. Since the first term in (15. 12ft is dominant near u = \/fkf, we 
arrive, in view of (I5.13I) , at the asymptotic linear growth law: 

a(t) ~ \JH 0 1, too. (5.16) 

It is interesting, and perhaps surprising, to note that makes no appearance here. 
Another interesting example that fits our study is the Friedmann equation 

H 2 = — — (pd + Pb + p 7 ), (5-17) 

considered in [22| in the context of the so-called geometrothermodynamics (GTD) [23]. 
where pd is the GTD dark matter density given in terms of the scale factor a by 

Pd = Pdo ( Aa ~ 3(Q - /J) + (1 - ^)a- 3(1+Q) ) ^ , d>0,a,^6l, (5.18) 

and pi, = Pbo a ~ 3 and p 7 = p 10 a~ 4 are the baryonic and relativistic matter densities. 
Thus, when some terms in these densities are absent so that the integration of (15. 1 7|) is 
of a binomial type, the Chebyshev theorem becomes applicable. Here we discuss such a 
situation. 

For simplicity, we assume the baryonic and relativistic matter densities in (15. 17[) are 
absent. Hence we arrive at the integration 


±2 



t = 


— 


J ((1 _ .4) + ^a 3 * 1+ '») do. 


(5.19) 


Thus, in view of the Chebyshev theorem, we see that, when a,/3 are rational numbers, 
la,p is elementary if and only if one of the following three quantities 


/3 — a 


2(1 + ex) 2(1 + (3) 2(1 + ex) 2(1 + /3) 2(1 + cx)(l + (3) 

is an integer. For example, when a — (3 — — |, we have 


(5.20) 


/_3 _3 = .Aina + -(1 — *4.)ez 2 + C, 


(5.21) 


where C is an integration constant. Hence we see that the big bang solution, a(0) = 0, 
is not allowed, and when A < 1 the scale factor obeys a power law growth pattern and 
when A = 1 an exponential pattern. 


5.2 A Two-Fluid Model 

Consider the equation of state of a two-fluid model 

2 R A 2 -R 2 . . 

Pm — — H- 4 —, (5.22) 

a 6 a 4 
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where A, R > 0 are constants and both cases A > R and A < R are of cosmological 
interest. If n = 3, the first term corresponds to pressure free matter and the second to 
radiation. We are interested in a closed universe with vanishing cosmological constant. 
In that case the Friedmann equation reads 


l6irG n ( A 2 — R 2 + 2 Ra — B n a 2 ) 
n(n — 1) a 2 


(5.23) 


where 


B r 


n(n — 1) 

167rG n 


(5.24) 


The necessary condition for (15.231) to be meaningful, that is, the right-hand side of (15.231) 
is positive for some values of a, is 


r>2 

R 2 < A 2 + 

BA n 


which guarantees that the quadratic function 


q(a) = A 2 — R 2 + 2 Ra — B n a 2 


(5.25) 


(5.26) 


has two real roots, which are 


a± = 


= 


B r 


—) +-(A 2 -R 2 ). 
B n Br 1 


(5.27) 


If A > R , then a + > 0 and a_ < 0; if A < R, then a + > a_ > 0. In either case, we 
may rewrite (j5.23j) as 


( _ ( a a -)( a + a ) 

«/ 9 > 

a 2 

where a stays in the interval 0 < a < a + if A > R] in a_ 
formally leads to its integration 


(5.28) 


< a < a + if A < R, which 


I = I a ((a — a_)(a + — a)) 2 da = ±f. 
Setting u = a — a_, we have 


R 


Br 


R 


— + —(A 2 — R 2 ) —— = —a- < u < a + — a_ 


B, 


B, 


2 ' n ir) +4 -(a 2 -&), 

BAn. / BArt. 


and we can rewrite the integral / on the left-hand side of (j5.29[) as 

/ = / a 5 ((a + — a_) — w)~ 2 du + a_ / u~ 5 ((a + — a_) — u)~ 2 du = I\ + Jo. 


(5.29) 


(5.30) 


(5.31) 
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It is easily seen that the Chebyshev theorem implies that both I\ and I 2 are integrable to 
yield elementary functions. Carrying out the integrations, we obtain 


/ = -(a+ + a_) arctan ( —— i = ] — ^yJa^-afi/l--aZ) + C, (5.32) 

2 \ 2y/(a+ — a)(a — a_) / 

where C is an integration constant. If A > R, the big bang solution, a(0) = 0, is allowed; 
if A < /?,, the initial value of the scale factor has to satisfy a(0) > a_. These fix the 
constant C. 

Below we consider the big bang solution with the assumption A > R. Thus C in 
(15.32)1 is given as 

C = Co = \J a_|_(—a_) H—(a+ + a_) arctan [ — + = ) . (5.33) 

2 V 2 v a +(- a -)/ 

It is clear that the solution increases monotonically until it hits a = a + when t becomes 

t — t s — Co + —(a+ + ci—). (5.34) 

From (15.28)1 we find a(t s ) = 0. Furthermore, differentiating (15.28(1 and taking the limit 
t —* tj, we have 

a(i.) = - (a+ 2 ~ a-) < 0. (5.35) 

This indicates that the solution starts to decrease beyond t s . In fact, if we use a(t) to 
denote the solution obtained above in the time interval 0 < t < t s , the solution of (j5.28)l 
that continues to evolve beyond t s is given by 

a(2t s — t), t s < t < 2 t s , (5.36) 

as can be verihed directly. At the time 2 t s , the solution vanishes again and then repeats 
the same cycle. That is, we have obtained a periodic solution which starts from a = 0 
at t — 0, climbs to its maximum height a + at half period ~t p = t s , and then starts to 
descend steadily until a = 0 at a full period t p = 2 t s . The quantity t s may be called 
time of maximum expansion. In fact t s is a moment of time symmetry since the solution 
curve of the scale factor against time is clearly symmetric about its time of maximum 
expansion. Typically such solutions expand from a “Big-Bang” and recollapse to a “Big 
Crunch”. In the case of a negative cosmological constant, the Big-Bang is a coordinate 
singularity and the solution may be extended to a periodic solution. I 11 the case of a 
barotropic equation of state, both the Big Bang and the Big Crunch are singular there is 
no unambiguous periodic continuation. 

For the record, we can explicitly rewrite in view of (15 .2 711 . (15.33)1 , and (j 5 ,34ft the period 
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t p just obtained as follows: 


7T 


t p — 2Cq + — (a++ a_) 



nGJA 2 - R 2 ) 2n 2 G n R AnG n R 

-7-77— + “7-77 + —( -77 arctan 

n[n — 1) n[n — 1) n[n — 1) 


ARV^G ri 


\Jn(n — 1)(A 2 — R 2 ) 


(5.37) 


Thus we have arrived at a comprehensive and accurate understanding of the evolution 
of the scale factor in the two-fluid closed-universe model with respect to cosmic time. 

The problem of this section motivates us to study the equation of state for the situation 
that the mass density and scale factor are related simply by 


Pm 


A B 

-1-■ 

a r a s 


r 7 ^ s, 


(5.38) 


where A, B , r, s are constants. Using (I5.38jl in the conservation law (I2.7jl which may be 
rewritten in the form 


d ® clp m 

m Pm i •> 

n da 

(5.39) 

we find 

a r \n ) a s \n ) 

(5.40) 

We are interested in the case when P m is simply a power function of a. 
loss of generality, we may assume 

Thus, without 

r — n, 

(5.41) 

which leads us to arrive at 


a s P m = B ^ - l) , s^n, 

(5.42) 


and the equation of state 


Pm 


( \ — 

TlPm \ s 

B{s-n) J 


(5.43) 


5.3 A Trinomial Friedmann Equation 

We consider the trinomial Friedmann equation with k — 1 and a nonzero cosmological 
constant A of the form [24j 


167rG n 2A 1 

n(n — l)^ m n(n — 1) a 2 ’ 


(5.44) 


where the mass density p m is taken to obey the specific relation 




















in which C\ is a constant, so that an integration of the equation leads to 


± t — I n — / -1 + 


2A o WnGnCi _ (n _ x) \ 5 


- a 2 + 


n(n — 1) n(n — 1) 


da. 


(5.46) 


For n — 3 and A > 0 (say), we have 


2 A 4 8 ttG 3 Ci\ 2 
a —a + —a H--- | da 




(setting u — ar) = -\/ — / ( tr —-w + 


1 / 3 


2 V A 


A 


87^361^ 2 

A~ 


d u, 


(5.47) 


whose integration is studied in [24] using the Weierstrass elliptic functions. Here we 
attempt a direct integration in view of Chebyshev’s theorem and in terms of elementary 
functions. For this purpose, we assume 


32nG 3 C 1 A < 9, (5.48) 

(the critical case when 32tiG 3 Ci = 9 is trivial) which allows us to rewrite (I5.47j) as 


h= yi [ - ui ^ u - u 2 dw ’ 


where 


Ul ’ 2 2A ± 2 V [ A 


327tG 3 Ci 

A 


(5.49) 


(5.50) 


are ensured to be real and u\ > u 2 in view of (I5.48j) . Setting v = u — u\ (say), we get 


/s 2 V A 


'2 (v + (u 1 -u 2 )) 2 dv, 


(5.51) 


which is integrable by the Chebyshev theorem. Hence, carrying out the integration (15.5 ip . 
we hnd 


h = 



-i(«i + u 2 ) + u + \Ju 2 - (ui + u 2 )u + UiU^j + C 


3 , , / 2 3 87rG 3 C! , 

’ 2A + “ + V “ _ A“ + A + C 


3 2 / 4 3 9 87r£r 3 Ci \ 

._ + 0 + J a - x a + X * + ^ 


(5.52) 


where C is an integration constant. 
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It is readily checked that the condition (15.48p does not permit a big bang solution with 
a(0) = 0. In fact, to stay in the real regime, it suffices to impose the condition 


a(0) > y/u[ = 



87rG 3 Ci 

A 


1 

2 


(5.53) 


for the initial value of the scale factor. Hence the expanding-universe solution following 
from (15.52p is given by the formula 


a 2 + 


3 2 , 8 nG 3 C 

T a 3-r 

A A 


I = 4. + C °e 2 ^, 


(5.54) 


where Co > 0 is a constant depending on a(0). Consequently, we have the following 
exponential growth pattern for the asymptotic form of the solution: 


a(t) ~ 



t —y oo, 


(5.55) 


which spells out the role of the cosmological constant A clearly. 

On the other hand, it can be examined that /„, when n > 3, does not allow a reduction 
into an integration involving a binomial differential so that the Chebyshev theorem is 
applicable. 


5.4 Friedmann’s Equation in a Chern—Simons Modified Gravity 
Theory 

Following Jackiw and Pi [25], the Chern-Simons modified gravity theory in four spacetime 
dimensions is based on the extended Hilbert-Einstein action 

s = i6b/{ y=5fl+ i‘ M } dl1 ’ (5 - 56) 

where R is the Ricci scalar, *RR the Pontryagin term given by 

•RR = -rr^Ri^, (5.57) 


with Rafj, v the Riemann tensor and *R^ V = i e^ vtl ' v 'R^ the dual Riemann tensor, 0 a 
prescribed quantity, and the cosmological constant is taken to be zero. Recall that the 
Pontryagin term *RR may be represented as a total divergence [231 26, [27]: 


\*RR = V P A^, 

where the 4-current reads 





(5.58) 


(5.59) 


30 












in terms of the Christoffel symbols, resembling a Chern-Simons invariant [ 28l [29] . 

In [30], a general Chern-Simons modified gravity model defined by the action 

S = ^ j | V=jR + J&-RR- i^ea„e] dV + s mat , (s.eo) 

is considered, where 0 is now a dynamical variable [31], of the flavor of the Bekenstein 
model [[32], [33, 34], l a coupling parameter, and 5]^ denotes an additional matter action 
contribution. In the Friedmann-Robertson-Walker metric limit governing a homogeneous 
and isotropic universe, the scalar held 0 is shown 1301 to obey the dynamics 

0 + 3HQ = 0, (5.61) 

which gives rise to the first integral 

0 = Ca~ 3 , (5.62) 


where C is a constant, and the energy conservation law reads 


^ (p + Y a ~ 6 ) + 3H ( p + P + ° 2a ~^ = °> 

which is seen to simplify into the usual one 

p + 3 H (p + P) = 0, 


(5.63) 


(5.64) 


independent of the Chern-Simons dynamical constant C since H — On the other hand, 
the hat-space (k = 0) Friedmann equation is 


H 2 = 


8 v tG 


C 2 

p + Y a 


-6 


which relies on C. Thus, with the linear equation of state 

P = (7 - 1 )P, 


(5.65) 


(5.66) 


we read oh from (I5.64p the relation p = p 0 a 37 . Consequently we arrive at the modified 
Friedmann equation 


87 tG 


—37 ^ C 2 -6^1 _ SvrG 

p ° a + T a = ~^~ PeS ' 


(5.67) 


where p e Q is the effective energy density given by 


C 2 


Peff deff ( a ) Po a ~h cl , 


(5.68) 


resembling another two-huid model. 
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Integrating (15.67H . we find 


I = j a 2 37 + da = ±2yj^^-t. (5.69) 

In view of the Chebyshev theorem, we see that, when 7 is rational, the integral I is 
elementary if and only if 7 assumes the values 

7 = 2-4 N — ±1, ±2,..., (5.70) 

7 = JV = 0,±1,±2,.... (5.71) 

As some illustrations, we first choose N among the sequence (15.701) . say N = 1 or 
7 = 1 (dust), to get 

1 = 37 + t ) + Cl ' (5 ' 72) 

where C\ is an integration constant. Hence the big-bang initial condition a(0) = 0 gives 
us the expanding solution 

a(t) = (QnGp 0 t 2 + 2V^G\C\t^ 3 , t > 0. (5.73) 


Similarly we choose 7 among the second sequence (15.711) . say = 1 or 7 = | (radiation), 
to get 


1 - 27 1 Po “ + ~ 


C 2 \* C 2 , ( 


V 


C 2 


In + I Po a + ) I + Ci, 


(5.74) 


where Ci is an integration constant so that for the big-bang solution it is fixed to be 

(5.75) 


4po 


From (I5.74p . we obtain the following asymptotic behavior for expanding solutions: 


a(t) ~ 2 


2nGp 0 \ 4 1 


3 J 


t 2 , t —>■ 00 . 


(5.76) 


So in both cases the solutions increase as power functions. 

Thus, we have seen that, in the multiple fluid situations studied above, power-law 
growth, periodic oscillations, and exponential expansion may all occur. It is interesting 
to note that a very wide range of evolutionary behavior can be achieved in the context of 
appropriate single-fluid models [35j as well. See also [36]. 
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6 The Reduced Temperature Equation 

In this section we study the Friedmann equation in terms of the reduced temperature and 
conformal time. We will see that, in such a setting, although the Chebyshev theorem is no 
longer applicable, new insight regarding integrability may be acquired which complements 
what may be obtained in terms of the scale factor and cosmic time. 


6.1 Governing Equation in Conformal Time 

We follow [57] to consider the normalized Friedmann equation 

1 /da \ 2 C r C m k A 
a 2 \ dt) a 4 a 3 a 2 3 ’ 


( 6 . 1 ) 


in 3 dimensions, where C r and C m are radiation and matter parameters, respectively. In 
the zero cosmological constant limit, A = 0, the equation can be integrated in view of the 
discussion in Section [5] In this section, we study what happens when we switch to the 
conformal time 77 : 

dt = adr], ( 6 . 2 ) 

and reformulate the Friedmann equation in the reduced temperature 


r = 


Vk = 


as m 


a/A 7 c a v L 3C r 

. In terms of these new variables, the Friedmann equation (16.11) becomes 


(6.3) 


drY 


A 


in) = “ T “ + Y - fcr2 + 3 


(6.4) 


where 

A = —, a = C r A c , (6.5) 

C 

are dimensionless parameters called the reduced cosmological constant and the reduced 
radiation parameter. It is interesting to note [37] that the Hubble parameter now assumes 
the form 

H = In a = -a/Y^- ( 6 . 6 ) 

dt d7? 

When a = 0, the integration of (j6.4[) may be carried out in terms of elliptic functions [37]. 
Besides, it is clear that, when A = 0 and a Y, Y 0, the question whether (16.4p can 
be integrated lies out of the reach of the Chebyshev theorem. Nevertheless, we show here 
that the latter case enjoys an exact integration as well as it is given in the cosmic time t 
and the scale factor a. 
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6.2 Integration of Equation 


We will only consider the nontrivial case, a/0. Thus we obtain from (16.411 the integral 
/ = ±77 where, formally, 


I 





1 ' 

a 

r H- 

3a 



1 

2 


dr, 


(6.7) 


which spells out another trivial case, 9 ak = — 1 . So we will only consider the nontrivial 
case, 9 ak ^ —1, below. Positivity of the quantity under the square root in (16.711 requires 
a > 0 or 

^- + k< 0 . ( 6 . 8 ) 

9a 

This observation splits our discussion into a few separate cases. 

« l ^ a + k > 0, then necessarily we have a > 0 in view of (I6.7H . Hence | + ak > 0. 
Now set 


Ti, 2 = 


-\±J\ + ak 


a 


Introduce the new variable v : 


((r - ri)(r - r 2 )) 2 = v{t-ti). 


(6.9) 


( 6 . 10 ) 


Then r = ( r x v 2 — r 2 )/(r 2 — 1) and ((r — ri)(r — r 2 )) 2 = — (r 2 — Ti)v/{ v 2 — 1). Then (16.7[) 
becomes 


y/al = j' 


T (( T - t i)(t ~ T 2 )) 2 dr 


= -2 

If k = -1, we have r 2 < T\ < 0 and 


dr 


nr - r 2 


( 6 . 11 ) 


1 r 2 1/1 

TlT 2 = — = - - 

a Ti a \ 6 


1 

- a 

9 


( 6 . 12 ) 


Hence we obtain 


/ = In 


3 ^/av — (l + y/l — 9a) 


+ C, 


(6.13) 


3^/a v + (l + -\A — 9a) 

where and in sequel C denotes an integration constant. If k — 1, then r 2 < 0 < T\ and 

7 1 = ~k(l + ^lr^ ■ (6 - 14) 
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So we obtain 


/ = —2 arctan 


3y/av 


+ C. 


(6.15) 


1 + y/1 + 9ct 

(ii) K £ + k < 0, then both subcases, a > 0 and a < 0, are allowed. Hence, for 

a > 0, we have k = —1. So, taking r + = u, u = \J—^ (—1 + <^) tan#, and tan | = v, 

consecutively, we obtain 


I = 


a 


T + 


3a 


1 

a 


-1 + 


9a 


dr 


= 6 yfa 


dt> 


= In 


V — V\ 


v - v 2 


+ Qa \J~a ( _1 + 5h) V ~ 1 

+ G, 


where 


v 12 = 3 y/a 


9a 


± 1 


(6.16) 


(6.17) 


are the roots of the denominator of the integrand of the second integral in (I6.16j) . which 
further renders (16.16j) into a more explicit form 


/ = In 


v + \/9 a — 1 — 3i/a 


+ C. 


(6.18) 


v \/9a — 1 “I - 3-y/a 

Similarly, for a < 0, both k — — 1 and k = 1 are allowed. So, taking T + ^ = u, 


u = 


\fa ft + sin#, and tan | = v, consecutively, we obtain 

/ = 6y/—a j - dV 

J 7 


(6.19) 


-6aVs i k + <k) v + 1 

Note that the discriminant of the denominator of the integrand of (j6.19[) is 

A = 36 ak. (6.20) 

In the nontrivial situation, k ^ 0, the condition a < 0 leads to A > 0 when k — — 1 and 
A < 0 when k — 1. Thus, we arrive at 


In 


v+\/l—9a— 3V —a 
j _ J v+Vl~9a+3V~a 

2 arctan v+ ^^ 1 + C , 


+ C, k — —1, 
k = 1. 


( 6 . 21 ) 


It is interesting to notice the similarities between (I6.13[) . (16 .1 8f) . and the first formula 
in (16 .21 [) . and between (16.151) and the second formula in (16.211) . 
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7 Friedmann Type Equations in Other Contexts 


It is worth mentioning that the Friedmann type equations are also encountered in numer¬ 
ous other physical situations out of the remit of relativistic cosmology. Below we list a 
few examples. By this means one may hope to extend the range of the applications of 
the Chebyshev theorem but more generally to gain insight into the underlying geometry 
behind the Friedmann equation and to spot useful ways of representing the solutions. 
One particular feature that emerges is that many of the solutions may conveniently be 
represented as roulettes. This will be the subject of a future paper. For the present we 
content ourselves with pointing out some correspondences. 

We begin by noting that if k — 1 we have encountered the Friedmann equation in 
three different but equivalent forms 


1 /da 
a 2 \ d£ 

ey 


/ dr 
\di] 


\ 2 1 

+ = PW, 

7.1 

j a- 

2 

+ CL? — 0 ? p(^) i 

(7.2) 

+ r 2 = p(i), 

(7.3) 


where dr/ = a and r = 1 


Note in this section we will adopt units in which 1 f 7rG ' 1 T \ = 1. 

^ n[n— 1 ) 


7.1 Refraction in a Horizontally Stratified Medium 


Snell’s law for the trajectory of a light ray y 
reads 



y(x) in a horizontally stratified medium 


K 2 y 2 ’ 



where K > 0 is a constant and n = n(y) is the refractive index. Equation (17.4p is clearly 
of the k = 1 Friedmann equation form with y(x) corresponding to the scale factor a(t). 
A direct integration gives us 


/ 



±x. 


(7.5) 


As an example, consider 

n(y) = Ay a , (7.6) 

where A, a are constants, d / 0. This corresponds to a barotropic equation of state with 
3y = 2(1 - a). 

Applying the Chebyshev theorem, we see that all the integrable cases for a are: 

« = ^, N — ±1, ±2,.... (7.7) 
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In particular, dark energy, radiation and dust correspond to a = 1, —1, — respectively. 
Explicitly we have 


/ 


A 2 


arctan 


A 2 

KS 


A 2 

K*y 


1 + C, 





OL = — — 
uc 2 , 

a = — 1 , 
a = 1 , 


(7.8) 


where C is an integration constant. 

Note that if a > 0 the refractive index increases with height, the rays are bent upwards, 
and this gives rise to a typical mirage effect. For example in the case of a = 1, i.e. dark 
energy, the rays are catenary curves. By contrast if a < 0, the refractive index decreases 
with height and the rays are bent downwards. If a = —1, i.e. for radiation, the rays are 
semi-circles and for a = — i.e. for dust, the rays are cycloids. In fact all three curves 
are roulettes. The case of dust, a = —| is equivalent to the brachistochone problem. 

Another interesting case is when the index n in (17.41) assumes the form 


/ \ Ay 

n{y) = , 

i + y 

(7.9) 

A 

(7.10) 

„ = ->!. 


where A is a constant satisfying 


The refractive index increases upwards and the rays are bent upwards. Hence the integral 
/ in (17.51) becomes 


/ = 


l + y 


vV - 1 )y 2 -2y-\ 


d y 




O- 


( 1 + y) 


(7 — 1 


y + 


<7 + 1 


d y, 


y > 


a — 1 


or y < 


C7+1 ’ 


(7.11) 


which in view of the Chebyshev theorem may be integrated similarly as done in Section 
6 and is omitted. Here, we simply write down the result: 


I = 


a 2 — 1 


V(a 2 - 1 )y 2 -2y-l + 


(7 


(u 2 -l)i 


In 


y 


a 2 -1 


1 +Jy°- 2v + 1 


a 2 — 1 


where C is an integration constant. 

A difficult case is when the index n is given by 


n(y) = 


C 0 y 2 


Co > 0. 


(1 + y 2 r 

This problem lies beyond the reach of the Chebyshev theorem. 


+ C, 

(7.12) 

(7.13) 
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7.2 Soap Films and Glaciated Valleys 

Soap films are modeled by surfaces which extremize surface area for fixed volume enclosed. 
This is equivalent to their being of constant mean curvature. Thus if the soap film is a 
surface of revolution obtained by revolving the curve y = y(x) about the x-axis we have 


5 


(Wl + fi/'j 2 - Ay 2 ) dz = 0, 


(7.14) 


where A is a Lagrange multiplier enforcing the constant volume constraint. 
Lagrange equation is 


+ 


(1 + y'T 2 a + y'T 2 


= 2A 


The Euler- 
(7.15) 


which is the statement that the surface has constant mean curvature since the principal 
curvatures are 


1 y" 1 _ 1 

Ri (1 + y ,2 )i’ R r i y(l + y /2 )5 


(7.16) 


Noether’s theorem gives a first integral 


-+ A y 2 = constant = c, 

\/i + y' 2 

which may be written in the Friedmann form 

\y ) y 2 ( c-\y 2 ) 2 

One has the correspondence 

1 

a++y, t^x,, pH--— 

(c- Ay-) 2 


(7.17) 


(7.18) 


(7.19) 


Note that in n spatial dimensions, when the 50(2) invariance is replaced by the SO(n— 1) 
invariance, Ay 2 must be replaced by y” _1 (see [45]). It was discovered by Delaunay Til '-Itl 
that we get in this way the locus of the focus of a conic section rolling without slipping 
on a line. There are three cases, the parabolic catenary, the elliptic catenary, and the 
hyperbolic catenary. The surfaces of revolution are called the catenoid, the unduloid and 
the nodoid 

If n = 3 and A = 0 we have a minimal surface of revolution which is well known to 
be a catenary. This corresponds to dark energy. If n = 3 and A 7 ^ 0 the solution may be 
given in terms of elliptic functions. If n > 3 we require in general hypo-elliptic functions 


The case c = 0 is rather pathological. Superficially, As a soap film it corresponds to 
a cylinder with y' = 0. This is a perfectly valid solution as a soap film . However its 
interpretation as a solution of the Friedmann equation it would correspond to a static 
closed universe supported by radiation with p = > 0. (This is a well treated case in 
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literature [Hj. See also [I] for a systematic treatment.) However a = 0 does not satisfy 
the Raychaudhuri equation 

~ = ~(P+3P) (7.20) 

a 3 

Alternatively if one sets H = 0 in equations (2,5) and (2,6) one obtains a contradiction 
unless |p + P = 0. 

In fact it is easy to show that the solution of the k — 1 Friedmann equation for 
radiation in terms of conformal time is given by 


a(v) = ^ cos r^, 


t 


1 

— sin 77 
A 


(7.21) 


which corresponds to a semi-circle in (x, y) space and hence a sphere when revolved around 
the x axis. 

The example above was based on a variational principle and may be further general¬ 
ized. The Euler-Lagrange equations of the variational problem 



<5 J L(y,y') dx = 0 , 

(7.22) 

where 

L = f(y)\/l + {y ') 2 -g(y) 

(7.23) 

admit the first integral 

f 

+ W 9 c ’ 

(7.24) 

whence 

( S V + 1 = /2 , 

y y 2 y 2 (c-y ) 2 

(7.25) 


An interesting example is provided by a theory of the shape of glaciated valleys [45] 
for which f — y and g = Ay, where the constant A is a Lagrange multiplier. Thus we have 


*d +1 = 

y 


y 2 (Ay - cf 


(7.26) 


That is, 


(v'Y = 


y 


- 1. 


(Ay - c ) 2 

First set u — Ay — c, assuming A 7 ^ 0 and c 7 ^ 0, otherwise it is trivial. Then we have 


(7.27) 


, 2 (1 — A 2 )u 2 + 2 cu + c 2 


(7.28) 


Thus we may assume A 7 ^ ±1 otherwise it is trivial. With A 7 ^ ±1, we use the transfor¬ 
mation 


w = (A 2 — 1)- — A 

c 


(7.29) 
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for comparison with the result in [35] to recast the equation (17. 2 7ft into 


a (A 2 - l) 3 (1 - w 2 ) 
c 2 (A w + l) 2 ’ 


(7.30) 


which can be integrated in view of the Chebyshev theorem. Carrying out the integration, 
we have 


± 


(A 2 -I)i 


-x 


c 


= — A\/l — w 2 + arcsinw + C, w 2 < 1 if A 2 > 1, 


(7.31) 


as given in [ [5j. and 

> 2 ^ 

I c 


± 


(1-A 2 )§ ,- ,- 

-—:- x = XV u> 2 — 1 + In w + Vw 2 — 1 


+ C, w 2 > 1 if A 2 < 1, (7.32) 


where C is an integration constant. 

The case c = 0 deserves special treatment. As a Friedmann equation it corresponds 
to a barotropic equation of state with p > 0 but p + 3 P = 0. It thus permits a solution 
for which a = 0. In general the solution is linear in conformal time. This is consistent 
with the equation of state just violates the strong energy condition. 


7.3 Catenary of Equal Strength 

This curve satisfies ([48] p. 305 see also [50] p. 192) 


y 


(i + (y') 2 ) c ’ 


(7.33) 


with first integral 


whence 


x 


arctany' = —b A, 
c 


y = B — c In 


cos ( - c + A 


(7.34) 


(7.35) 


Inverting (j7.34j) for y' to get ( y') 2 as a function of cos 2 (- + A) and using (17.35ft lead to 


i fy'Y i f2(y-B) 
72+ - = T7 exp 


y \y 


y 


(7.36) 


which is of the Friedmann form under the correspondence (x,y) O (t,a). Viewed as a 
closed Friedmann model it is time symmetric and exhibits a “Big Rip” [IS], i.e a blow up 
of the scale factor, at finite time in the future and the time reverse of a big rip at finite 
time in the past. However if one evaluates the pressure P the resulting equation of state 
does not seem to be very simple nor to have a simple physical interpretation. 

The equation (I7.36jl is obviously beyond the reach of the Chebyshev theorem. 
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Identifying the “matter density” p as 


1 2y 

P = — e c 

y 

and applying the “energy conservation law” 

p' + Z(p + P) V - 

y 


o, 


(7.37) 


(7.38) 


we deduce the equation of state relating the “pressure” to “matter density” as follows 

P = -\p- 7 Tpy = -\p- 7 rpf(p), (7-39) 

3 3c 3 3c 

where / is the inverse function defined by (I7.37j) if it exists. This equation is complicated. 
As a by-product, we have p + P = |p(l — |) = |p(l — ^—) whose sign is undetermined. 
However p + 3P = — — and therefore if y > 0 and c < 0 we have p+P > 0 and p+3P > 0 
and both the dominant and the strong energy conditions hold. On the other hand, if c > 0 
then the strong energy condition is violated. 


7.4 The Elastica of Bernoulli and the Capillary Curve 

Both of these curves satisfy (see |50] p.275) 


± 


y 


(1 + (; y 0 2 * ) 


with the first integral 


Rearranging (17.4111 we have 


r 


T 


y 

\ b 2 ’ 

(7.40) 

y 2 

2 b 2 c ' 

(7.41) 

“Hydrostatic Arch”. 



V 2 (c + & 


2 > 


(7.42) 


which is of the Friedmann form under the correspondence ( x,y ) yy (t,a). In general, the 
solution is given by an elliptic integral but if c = 0, ±1, then there are some elementary 
solutions (see below). 

The special case c = 0 has a linear equation of state corresponding to stiff matter with 
positive energy density and P = p H The Friedmann equation in terms of conformal time 
becomes 

+ 4a 4 = 166*, (7.43) 


2 The same curves arise in the case of the Mylar balloon [52]. If one interchanges x and y one obtains 

(17.4111 and the Mylar balloon corresponds to setting our integration constant c to zero. More generally, 

any equation of the form 


vh+b/T 


= f{x) may be converted to the Friedmann equation form in this way. 
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whence 


(7.44) 


a = a( 0 )\/cos( 2 ? 7 ) . 

Cosmic time is then given by an elliptic integral. 

If c = —1 and we take the upper sign in (I7.4ip we can find a complete solution for y 
and x in elementary terms. For a direct evaluation using the integrals see below. In order 
to better understand the global behavior of the solutions, we first recall the traditional 
approach (see e.g. [50]) and we will take contact with the discussion based on direct 
integration below. We define 0 to be y' = tan0, so 0 is the angle the curve makes with 
the rc-axis. the first integral (I7.4ip then becomes 


c =F cos 0 



(7.45) 


Setting c = — 1 and taking the upper sign in (17.451) we have 


y = 2 b 




(7.46) 


where we have chosen the positive sign for the square root and we shall ultimately be 
interested in 0 ranging from 0 to 2n. Now 


dx dx ds 
d 0 ds d 0 ’ 


(7.47) 


where ds = a /1 + (y 1 ) 2 dx along the curve, and using the fact that the curvature of a 
curve is given by 

y" d 0 


(i + (z/') 2 )§ ds ’ 


(7.48) 


and (I7.35P we find that 


da: , 1, cos -0 , 1, 

- = i —0-- = it — 0 

d0 2 sin t 2 


Sill 


9 • ^ 

-— 2 sm — 

t 2 


(7.49) 


Taking the upper sign in (I7.49P we see that x increases with increasing 0 from 0 = 0 to 
0 = 70 it then decreases with increasing 0 until 0 = ^ and thereafter increases with 
increasing 0 until 0 = 2n. Integrating (I7.49p we find cosmic time to be given by 


t — x — b In 



+ 2 b cos 


0 


(7.50) 


We see that as 0 runs from 0 to 2tt, x runs from — oo to +oo, and the curve, which is 
symmetric about the y-axis has a node on the y -axis and a loop running from the node 
and back thereafter continue to infinity. Between 0 = | and 0 = ^ the function y(0) 
is double valued. On the branch from 0 = 0 to 0 = n the function y(0) is monotonic 
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increasing from zero and on the branch from 0 = 7r to 0 = 27r it is monotonic decreasing 
to zero. The maximum value of y is 2b. At the node y = 2b sin where 0 n satisfies 

In ^tan + 2 cos = 0. (7-51) 

Considered as a closed Friedmann model conformal time is given by 


d^ = 


dx 

y 


dt 

a 



That is, with a choice of integration constant, 



(7.52) 


(7.53) 


If we take the branch that runs from t = — ootot = 0we obtain a cosmological model 
expanding from zero size in the infinite past and reaching a finite size in finite time at 
t = 0 at the node. The question then is what happens next? The continuation past the 
singular point is not obvious. If joined to the lower branch, then Hubble’s constant H = - 
would jump from a positive to an equal negative value at t — 0. However the Einstein 
equations would be violated and the jump in Hubble’s constant would correspond to a 
spacelike hypersurface on which some matter instantaneously appears from nowhere. 

A better alternative would be to follow the loop. The question is then what happens 
when one reaches -0 = |? Hubble’s constant is given by 


1 dy d0 _ 1 sin 2 | 

ydirdx cosf (1-2 sin 2 1)' 


(7.54) 


which diverges as 0 —)■ jj. This seems to be an example of a sudden singularity [15 lj. 
From f!7.42p we can read off the “matter density” to be 


P 



Hence with H 


— we have 
y 


P 


-2H 1 + 


(b 2 c + \y 2 ) 


P- 


(7.55) 


(7.56) 


Comparing this result with the conservation law (I7.38j) . we obtain 


:P- 


P = 


2 py 2 


3 ( b 2 c + \y 2 ) ’ 


(7.57) 
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which is positive for any c > 0. This is the strong energy condition, which implies the 
dominant energy condition p+P >0. If c < 0, the sign of p+P is generally undetermined. 

One might also consider the top portion of the loop. This gives a time-symmetric 
closed universe which starts with a big bang at finite size y = V2b, with an infinite 
Hubble constant and reaches a time of maximum expansion at a finite time thereafter 
after which it recollapses to a big crunch. 

In both of the two cases just discussed, the singularity is clear from (I7.42[i since if 
c = — 1 the energy density diverges as y —» y/2b as . This is reminiscent of the 

behavior of an equation of state with P = — |p . 

Similar possibilities are possible for solutions for which c ^ -1. In particular if c is 
negative similar singularities can arise. If c > 0 the behavior is more like the c = 0 case. 

In general, setting 

«=|^ + c (7.58) 

in (17.42p leads formally to the integration 


\/2 f u du _ 

N J \J{u c) (1 u 2 ) 


(7.59) 


which lies beyond the reach of the Chebyshev theorem unless c = 0, ±1. So it is worthwhile 
to boil down the range of c in which the problem is relevant. For this purpose, we rearrange 
(17.42}) to get 

(-^>4- cl -3( i+c+ ®' (7 - 60) 

So we see that the value c > 1 is prohibited and c = 1 gives rise to the trivial solution, 
y = 0. For the remaining nontrivial integrable case, c = —1, we have 


/_ i = \f 2 arctanh 


= v^arctanliA/l 


u 


2y/l -u + C 


y 

4 b 2 


2V2\ 1 


y 

4 b 2 


+ 


(7.61) 


where C is an integration constant. To make contact with our previous discussion we set 


and use the identity 



arctanh cos — = In cot —. 

2 2 


(7.62) 


(7.63) 
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7.5 Central Orbits 


Consider a particle moving in a plane with polar coordinates r, 9. If v is the speed and 
V(r) the potential energy per unit mass then by energy conservation, 

i(f 2 + r 2 0 2 ) + V(r) = S, (7.64) 

where £ is the energy per unit mass. Angular momentum conservation gives 


r 2 9 = h, 

where h is the angular momentum per unit mass. Thus 

1 /da\ 2 1 2 , r/ , . 

^(d^j + 72 “ 

Now setting u = - we obtain 



(7.65) 


(7.66) 


(7.67) 


We see that this equation is of the same form as the reduced temperature version of 
the Friedmann equation (16.4ft with k — 1, under the correspondence (u,9) -(-?• (r,r}) and 
bearing in mind that while 6 is a periodic coordinate period 27T, conformal time r/ is not 
periodic. 

Central orbit problems arise from considering geodesics in spherically symmetric met¬ 
rics. In this connection it is perhaps of interest that the equation (16.41) for the reduced 
temperature is of the same form as that obtained for null geodesics in the Reissner- 
Nordstrom metric (see equation (51) in [38]) except the sign of the cosmological term is 
negative for real values of the electric charge. 


7.6 Spherical Symmetric Lenses 

A closely related example to the previous one is provided by applying Fermat’s principle 
to a lens whose refractive index n(r) depends on the radial coordinate r. The motion may 
be assumed to be in a plane through the origin and we have 


n 2 (r) ( r z + r z 9 l ) = 1, 


(7.68) 


and 


n 2 r 2 0 — h — nr sin i, 


(7.69) 


where h is a constant and i is the inclination between the ray and the radial direction. 
Thus 


dr 

r 2 dd 


1 7T 

r 2 Jj2 


(7.70) 


45 







and so if u = - we have 


2 



+ iT 



(7.71) 


Note that while the lens problem is defined in the plane and hence 9 is periodic of 
period 2i r, no such restriction is placed on t. Thus the cosmological problem is defined 
on the universal covering space of the punctured plane. 

Some examples are as follows: 


Maxwell’s fish-eye lens is given by 


n(r) = 


n 0 


and thus 


1 + r 2 


Tin 


P = 


(7.72) 


(7.73) 


h 2 (l + a 2 ) 2 

Comparing with (17.18jl and Delauny’s construction of soap films, it would appear 
that this case corresponds to the locus of the focus of an ellipse rolled without 
slipping along a line (i.e Ac < 0). Now the optical metric of Maxwell’s fish eye lens 
is that of a round 3-sphere pulled back to M 3 by stereographic projection and hence 
all rays are circles. If these circles enclose the origin and lie in a plane through the 
origin they presumably map to the elliptic catenaries on the covering space. 


• The Eaton lens has 


12 a „ 

n= \/ -1, 

r 


0 < r < a; n — 1, r > a. 


(7.74) 


Thus for r < a, the solutions of which are portions of ellipses 

- = 1 + e cos(</> — </> o) 
r 

with focus at the center r = 0, semi-latus rectum l, and eccentricity e, where 

h 2 


(7.75) 


l = 


e = \ 1 


b 2 
72 ' 


a v a 

The ends of the semi-major axes of an elllipse are situated at a radial distance 

l 


(7.76) 


1 — e 2 


(7.77) 


from the focus at r = 0. In our case one finds that this distance is a-independent 
of the impact parameter h. Because the refractive index is continuous across the 
boundary, rays entering the lens initially continue in the same direction are thus 
initially at one end of the semi-minor axis and initially moving perpendicularly to 
the semi-minor axis of their elliptical trajectory within the lens. It follows that they 
exit from the other end of the semi-minor axis in the opposite direction from that 
which enter the lens. Thus the Eaton lens behaves like a retroreflector: Any ray 
entering the lens exits in the opposite direction. 
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• As pointed out in [47J the Eaton lens is dual in the sense of Arnold-Bohlin to the 
Luneburg lens for which 



(7.78) 


r < a; n — 1, r > a. 


8 Conclusion 

In this paper we have discussed the explicit integration of Friedmann’s equations for equa¬ 
tions of state more complicated than that of a single component fluid of linear barotropic 
form. In particular we have made use of a theorem of Chebyshev which gives sufficient 
conditions for explicit integrability. However we also found examples, while which are 
outside the scope of Chebyshev’s result, nevertheless allow explicit integration. As well 
as the traditional form of the Friedmann equation is in terms of cosmic time, we have 
treated the case when the independent variable is conformal time. Less traditionally we 
have also considered the inverse of the scale factor (the so-called reduced temperature) as 
a function of conformal time. We have seen through a thorough study of the Chaplygin 
fluid model that integrable cases may often provide sharp insight into nonintegrable cases 
as well. In particular, based on some integrable examples, we have derived a universal 
formula for the asymptotic exponential growth rate of the scale factor in cosmic time, 
regardless whether the Friedmann equation falls in its integrable regime, and the formula 
reveals the coupled roles played by the conventional linear matter and nonconventional 
nonlinear matter, with regard to the presence of dark matter. Finally we have shown that 
all three forms of the Friedmann equation occur in areas of physics not directly connected 
with cosmology. This suggests that methods of dealing with the equation developed in 
those areas may transfer to the cosmological setting. In particular in a forthcoming paper 
we will illustrate this point by expressing its solutions as roulettes. 
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by National Natural Science Foundation of China under Grant No. 11471100. 
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